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We analyze the behavior of an ensemble of inertial particles in a one- dimensional smooth Gaussian 
velocity field, in the limit of large inertia, but considering a finite correlation time for the random 
field. The amplitude of the concentration fluctuations is characterized by slow decay at large inertia 
and a much larger correlation length than that of the random field. The fluctuation structure in 
velocity space is very different from predictions from short-time correlated random velocity fields, 
with only few particle pairs crossing at sufficiently small relative velocity to produce correlations. 
Concentration fiuctuations are associated with depletion of the relative velocity variance of colliding 
particles. 
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I. INTRODUCTION 

Starting from the work of Deutsch it has been 
known for some time now, that inertial particles in a ran- 
dom velocity field undergo clustering phenomena. These 
behaviors are observed in numerical simulation of turbu- 
lence as well (see e.g [1] and references therein) and are 
thought to give an important contribution to coalescence 
phenomena, e.g. in the process of rain formation [1, 
Both in the case of random fields and of real turbulence, 
spatial inhomogeneity of the statistics, though contribut- 
ing to particle segregation does not appear to be an 
essential factor. 

Random velocity fields with various statistical proper- 
ties, have been used to model inertial particles in turbu- 
lent fiows, both in the presence of gravity @, 01 and in 
its absence . Simplified models disregarding in part 
or in total the spatial structure of the velocity field have 
been introduced as well (see e.g. [13, El), to cope with 
the difficulty of the analytical treatment of the multi- 
particle statistics. 

In the simplest instance, an inertial particle in a tur- 
bulent flow is characterized by the relaxation time of 
its velocity relative to the fluid: the Stokes time ts = 
2/9 a^A/z/Q, where a is the particle radius, A is the ratio 
of the particle to fluid density (assumed large) and i/q is 
the kinematic viscosity of the fluid [12]. Experimental 
data [l3l | and numerical simulations [14l | both indicate 
that clustering is stronger for particles with rg of the or- 
der of the Kolmogorov time scale. This observation, and 
the fact that Stokes times for most atmospheric aerosols 
of interest lie in this range [S] have motivated substantial 
analytical effort in the small Stokes time regime, with 
smooth velocity fields mimicking turbulence at the Kol- 
mogorov scale P, [l^ . 

Some analysis has been carried on also for Stokes times 
in the turbulent inertial range [1, [igI . [l7j , but numerical 
simulations [18] indicate that concentration fiuctuations 
do not follow simple scaling rules. 

Due to the relative ease of analytical treatment, an- 



other regime which attracted attention is that of inertial 
particles in a smooth random velocity field, with correla- 
tion time te much shorter than the geometric scale r^/cru 
obtained from the correlation length and the ampli- 
tude (T„ of the velocity field fluctuations [l^, [13] ■ How- 
ever, this is a short-time correlated regime that is very 
different from the situation in realistic turbulent flows: 
TE(Ju/rv = 0(1). 

The analysis of random velocity field models has al- 
lowed to identify at least two clustering mechanisms, ex- 
pected to be present also in real turbulence. The first one, 
originally proposed in is preferential concentration 
of heavy (light) particles in the strain (vortical) regions 
of the flow. This effect has been recently observed also 
in the case of particles with very small inertia [H, [1^ . 
The second mechanism, already present in one dimension 
(ID) [2l|,[23|, is that of the particles catching one another 
in their motion, as they slip with respect to the ffuid. 
In the case of smooth velocity fields, this leads to the 
formation of caustics in the instantaneous concentration 
field, which, in turn, act as back-bones for the clustering 
process [1^. The clustering itself is a long-time process, 
which can be described in terms of the Lyapunov expo- 
nents of the particle pair dynamics (25l . [26j . 

Focusing on the second mechanism, it appears that 
caustics formation tends to become maximum for 
Ts/TE,TEUu/rv 1 [13, and that, at the same time, 
clustering becomes weaker. For larger rg, one expects 
that the particles be scattered by the velocity fluctuations 
they cross in their motion (28j as if undergoing Brownian 
diffusion, resulting in vanishing particle correlations and 
in clustering destruction. It is to be stressed that the in- 
terest of this limit is by no means academic, as particles 
in turbulent flows for which rg lies in the inertial range 
(or above) , will see smaller vortices precisely in this way. 

Purpose of this paper is to understand in detail how 
and under which conditions the uncorrelated limit de- 
scribed in HI] is achieved. It will appear that clustering 
destruction occurs in a very non trivial way, requiring 
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consideration, among the other things, of how particle 
correlations decay at scales comparable with or above 
Tv We are going to show that, provided we are away 
from the short correlation time regime te ^ cfu/r^, clus- 
tering at Ts 3> Te will not be dominated by the small- 
separation particle pair dynamics (small with respect to 
Tv) described in 21]. Particle pairs remaining close long 
enough for a Lyapunov exponent approach to be appro- 
priate, are still present, but their contribution to concen- 
tration fluctuations is negligible. 

The key mechanism for the destruction of concentra- 
tion fluctuations will appear to be that, at large ts/te 
and finite tecFu/tv, only particles with increasingly small 
relative velocities (but not small enough for a local the- 
ory on the lines of [111, [2^ to be valid) have a chance to 
be correlated. This is to be contrasted with the picture 
in (2?! , of correlation fluctuations disappearing as caus- 
tics occupy in the above limit larger and larger portions 
of space, and with the one in ^26.] of saturation to the 
space dimension of the particle distribution correlation 
dimension. Neither in [26l | nor in [27l |. however, was a 
quantitative prediction on clustering decrease provided. 

The present analysis will also allow to show that, for 
large ts/te and rather generic (not too large) values of 
<^uTE/rv, clustering is associated with smaller typical rel- 
ative inter-particle velocities, compared to what would be 
observed in the absence of correlations. Thus, there are 
circumstances under which, clustering may hinder rather 
than enhance coalescence phenomena. 

This paper is organized as follows. The main defini- 
tions and approximations will be presented in Section II, 
following by and large the notation of [l9l . [20I [2l|. In 
Section III, the pair particle dynamics will be analyzed 
in the large rg limit, identifying the relevant time and 
velocity scales for clustering. In Section IV, the relation 
between particle-particle velocity depletion and cluster- 
ing will be established, and heuristic estimates for the 
concentration correlation will be obtained, starting from 
the Fokker-Planck equation for the distribution of the 
relative particle velocity and coordinate. In Section V, 
a proof for the dominance of slowly approaching particle 
pairs in cluster generation will be provided. Section VI 
contains the conclusion. 



II. MODEL EQUATIONS 

Consider an ensemble of inertial particles transported 
by a ID zero mean, Gaussian random velocity field u{x, t) 
with correlation 



{u{x,t)u{Q,Q)) = cTlF{t/TE)9{x/r^), 



(1) 



where 5(0) = F{Q) — 1, the function g is assumed 
to be smooth and to have all necessary moments, and 
F{t)At = 1. The particles are immaterial, so that 
they can cross without interaction, and their velocity v 
is taken to obey the Stokes equation: 



Introduce the Stokes and Kubo numbers: 

S = ts/te, K = (TuTE/r^ 



(3) 



and choose units so that cr^ = T5 = 1; therefore: 

TE^S-\ r, = {KS)-\ (4) 

In the regime S" ^ 1, the particle velocity will change 
little on the lifetime of a fluid fluctuation and Eq. ([2|) 
could be approximated by a Langevin equation. This, 
independently of K , i.e. of the random field being short- 
time correlated or not. It is possible to substitute into 
Eq. dm: 



Tp{v) 



At 



(5) 



(6) 



with ^(i) white noise: (^(t)^(O)) = 6{t) and 
{u{x{t),t)u{x{{)),G)\v) 

the correlation time of the fluid velocity sampled by a 
particle moving at speed v. In the above formula, the 
averages are calculated along trajectories, with {.\v) in- 
dicating the condition on v. The average on the trajec- 
tory in the denominator, due precisely to the presence of 
concentration fluctuations, does not necessarily coincide 
with the space average {u^) = a^. 

From Eq. ([T]), we can put Tp ~ te, provided r^/ay S> 
Te, with the particle velocity variance. If Tp ^ te, we 
can estimate from Eqs. ([2]) and ([5]): 

al 0, S-\ (7) 

which is satisfied if ^ S. Following the same line 
of reasoning (see Appendix), it is possible to show that 
prescription problems in the definition of the white noise 
arise only at 0{K^/S). Notice that ^ S and S > 
1 imply X — K/S <^ 1, where x is precisely the force 
magnitude parameter introduced in 21|. This condition 
is easily satisfied unless one chooses to work in a frozen 
turbulence regime K ^ I. 

Turning to the relative motion of particle pairs, let us 
introduce difference variables 



V2 -Vl, 



X2 - Xi, 



where 1 and 2 label members of a particle pair, and indi- 
cate ui,2(i) = u(xi^2it),t). In the regime 3> 1, starting 
from Eq. ([2]), we can approximate the equation for the 
relative motion of particles with the Langevin equation: 



(8) 



v = Ts'[-v{t) + uix{t),t)]. 



(2) 



For S 3> K^, we have for the noise square amplitude 
B{r) = b'^{r), from Eq. (P): 

Bir) - IZiiMt) - u,mu2{0) - ^.i(0)]|^,r)dt 



c.4^-i[l-5(^Av)]. (9) 

The decoupling of the difference variables {v, r) from the 
center of mass variables 5(2^1 +2:2) and ^(fi -I- W2) de- 
scends basically from the condition Tp(v) te- 
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FIG. 1: Typical particle pair trajectories for e <^ 1 (a) and 
e > 1 (6). The exit time to the right is r ~ e^^^^rs. The 
small rectangles are space-time correlated regions for the fluid, 
of size r-v and duration te. The time rg quantifies the particle 
memory. The particle approach in (a) is governed by an 0(e) 
negative Lyapunov exponent |2l|. 



III. TWO-PARTICLE DYNAMICS 

In order to determine how concentration fluctuations 
are generated, we need to understand the particle pair 
dynamics at small separations. In the large S limit, we 
expect that the velocities of colliding particles be uncor- 
related, however, this condition turns out not to be au- 
tomatically satisfied. A condition for the particles to be 
weakly correlated, despite being close, is that the Stokes 
time be much longer than the time spent at |r| '^r^^ so 
that the velocity of the two particles will be the result of 
many contributions by uncorrelated regions of the fluid 
(see Fig. [T|). Under these conditions, particles will move 
ballistically at scale \r\ < r^, and they will cross at veloc- 
ity 1/ (jy = 5^^/^, after a time t ^ r^/ay ~ (i^^S*)"-^/^. 
The weak correlation condition is therefore r^/ay <^ ts, 
i.e.: 

e^K'^S:^!. (10) 

The parameter e is the same as in the theory in [iqI . [20| , 
who considered in detail the small e regime. For e <C 1, 
the time spent by a pair of particles at |r| < Tv before col- 
lision, would be longer than ts, and the particles would 
have memory only of that last portion of their history, 
when they were at |r| <C r^; therefore, v could not be 
treated, at the time of crossing, as a difference of uncor- 
related velocities [1^ . 

Notice that the conditions e <C 1 and S ^ 1 can 
be realized only if K S^^/"^ ^ 1, corresponding to 
a short-time correlated random velocity field. In this 
regime, e plays the role of an effective Stokes number 
(see in particular [1^), in which the place of te is taken 
by / [g'^te) , that is the diffusion time of a (non-inertial) 
tracer across a distance r^. Notice also that the condi- 
tion Ts ^ tetT-jI'^u, i-e. KS^ S* » 1, is not sufficient by 
itself to guarantee weak correlation, that in fact is not 
fulfilled in the range S~'^ < if < 5"^/^, where e < 1. 

Back to the large e range, although most particles 
will be weakly correlated, those approaching at speed 



\v\ ^ cr„ may stay at |r| <r^ sufficiently long to end up 
being strongly correlated. It is then important to under- 
stand whether an effect like clustering, which depends on 
pair-correlation, is due to the effect of the many weakly 
correlated particle pairs, or to the one of the few strongly 
correlated ones. 

We will answer this question in Sec. V. For the mo- 
ment, we identify the important time and velocity scales 
involved in the dynamics of slowly approaching particle 
pairs. We rescale variables: 

i=t/ty, i> = jy/jy^, r = r/r^, (11) 

where 

U = (2e)-l/^ V. - 2i/3e-i/6^„ (12) 

and Tv is given in Eq. ([?]). In the new variables, at 
short separations, the noise amplitude ^ can be Taylor 
expanded and Eq. ([5]) becomes: 

b = -{2e)-^l^v + fi; f^v. (13) 

Thus, iv and arc in some sense the time and velocity 
scales at the correlated scale r^. Notice that the condi- 
tion for Eq. (|13p to be meaningful as a stochastic differ- 
ential equation is ^ te, i.e. again x — K/ S ^ 1. 

To understand better the meaning of tv and fv, con- 
sider the separation of a pair of particles such that 
f(0) = and |i'(0)| ^ 1. This is equivalent to looking at 
the approach of a particle pair, with conditions imposed 
on the velocity at the moment of crossing, rather than at 
a previous initial time. 

In the initial phase, the separating particles will move 
at almost constant relative velocity, and f ~ The 

variation of v will be At> ~ i>(0) /q t'^(t')di' so that 
(Ai>2) ~ £'2(0)£3/3. The initial phase wiU last until 
Ai> ~ ^'(0), i.e. until i ~ 1, when f ^ £"(0). We then see 
that ~ 1, i.e. v — \s the minimum relative velocity 
to which the approximation of constant relative velocity 
at the correlated scale applies, and is the associated 
exit time: |7'(iv)| r^. 

The following phase will be determined by the linear 
nature of Eq. ^] for t < ts, the drift ~{2e)-'^^^i' 
can still be disregarded, and we expect v{t),r{t) ~ 
i>(0)exp(i). After a time of the order of — ln£>(0), the 
particles will have moved apart by f ~ 1, i.e. they 
will be moving out of a correlated region for u. In 
this final part of the separation process, Eq. will 
cease to be valid, with v obeying instead i) ~ For 
Ts 3> i 3> — lni>(0), V will behave like a Brownian mo- 
tion: v{t) ~ <^/^, while f will undergo a Richardson- like 
diffusion process: f (t) ~ P/"^. 

In conclusion, we have two groups of particle pairs. 

1. The great majority, for which, at the moment of 
crossing: v a^; their motion in a correlated re- 
gion can be approximated as ballistic and their per- 
manence time is T ~ e~^/'^Ts- 
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2. A smaller fraction, for which, at the moment of 
crossing: Iv] < v^, with ~ e~^^^ay; their dy- 
namics in a correlated region is determined by fluc- 
tuations of u{x,t) at scale < and their per- 
manence time grows logarithmically as u — » 0: 
t{i^) ^ e-^/Hn{iy,/\iy\). 

The fact that for e ;» 1, iv <C T5, has the consequence 
that only particles meeting at exponentially small rela- 
tive velocities 

~ cr^e""^/^ exp(— const. e^/"^T/r5), 

will have a chance to stay within a correlated region for a 
time T >Ts. An exponentially low level of molecular dif- 
fusion is thus sufficient to destroy these pairs. Even in the 
absence of molecular diffusion, however, the relative mo- 
tion of these pairs is characterized in the limit t/ts — > 00 
by a positive Lyapunov exponent 21| and no long-time 
contributions to clustering are present. The remaining 
particle pairs can be shown not to produce singular con- 
tributions to the two-particle probability density function 
(PDF) for 1/ [s^l- We can then estimate the fraction of 
particle pairs staying close for a time r, i.e. the frac- 
tion of particles crossing at speed Vt (provided of course 
T > ty), as 



^1/6 



exp(— const, e^^^t/ts)- 



This, again, is exponentially small for t > ts and large 
e, while the fraction of particle pairs that are merely not 
ballistic, is the much larger value Vv/a 

The situation in the small e regime is very different, as 
Eq. (|13p leads to instantaneous equilibration of the parti- 
cle velocity: (i^^lr) ~ CT^(r/r^)^ and the ballistic particle 
fraction associated with caustics formation, disappears 
[2l|,[H|33|. Hence, a-'^{iy'^\r = 0) ^ and the particles 
behave for e <IC 1 as a mono-disperse phase (although 
with velocity not locally equal to that of the fluid, as it 
happens instead in the small S limit [iBl)- At e ^ 0, all 
particle pairs are therefore co-moving. 



IV. THE PAIR DISTRIBUTION 

We have seen at the end of Section II that in the regime 
S ^ 1, S ^ K^, the difference variables z^, r decouple 
from the center of mass ones. This means that at equi- 
librium, the two-particle PDF will be in the form 



p(wi,2,a;i,2) = ^ 



V2 



where 17 is the length of the domain for 3- (Unless am- 
biguous, we do not use subscripts to identify PDF's re- 
ferring to different stochastic variables). Multiplying by 
iV^, with N the total number of particles in the domain 
VL and integrating over duidw2, we obtain the expression 
for the concentration correlation: 

(n(r)n(O)) = n'^n ( dvp{v,r) := fi^[l + /(r)], (14) 



where n = N/ Q, is the mean concentration. The quantity 
/(r) = ^p{r) — 1 gives the strength of the concentration 
fluctuations. 

The difficulty of the problem, compared with the e ^ 1 
limit, is the finite width in v of the PDF p{v,r). In the 
e <g; 1 limit, small r implied small z^, a long permanence 
at that particular r and the possibility of linearizing the 
Fokker-Planck equation in r along the lines of Ji, Ji^. 
In our case this is not possible, and we must consider the 
full Fokker-Planck equation, which, from Eqs. ([8]) and 
([9|), reads, at steady state: 



1 



vdrpiv.r) - d^[vp{v,r)] = -B{r)d^p{v,r) 



(15) 



(see |32| for another example of kinetic treatment of the 
particle pair statistics in turbulent flows). 

Taking moments of the Fokker-Planck equation (fT5|) 
allows to draw important conclusions on the connection 
between clustering and decrease of the relative velocity 
variance. The first two moments of Eq. ([TS]) read: 



(16) 



dr{{v\r)p{r))={), 
dr{{v\)p{r)) + {v\r)p{T)^Q, 



where {v'P\r)p{r) — ^ dvv^ p{v,r). We can impose bound- 
ary conditions to Eq. (jl6[) at r — > 00, where the solu- 
tion, corresponding to uncorrelated particles, is known. 
The first moment equation tells us that {v\r)p{r) is con- 
stant; but, at \r\ 3> correlations are absent; therefore 
p{r — > 00) = $7"^ and {v\oo) — 0; thus {v\r)p{r) — 
for all r. Substituting into the second one, we get 
dr{{v'^\r) p{r)) = 0, and, using again p{r — > 00) = fl^^ 



and {1/ \oo) 



2S^\ we get: 
r-^{,.^\r) = [l + f{r)]-\ 



(17) 



In other words, clustering in ID is necessarily associated 
with velocity variance decrease. Notice that this result 
rests only on the applicability of the Fokker-Planck equa- 
tion which requires 5 » l,x = K/S 1, but is 
otherwise independent of e. 

An idea of how /(r) goes to zero as e — > 00 can be ob- 
tained treating the function g in Eq. ([S]) as a perturba- 
tion. This is motivated by the observation (see previous 
section) that most particles spend very little time at dis- 
tances |r| < Tv, and by the assumption (verified in {34|) 
that the remaining particles produce only a correction to 
p{r). To lowest order, the particles are therefore inde- 
pendent random walkers with velocity correlation time 
Ts- A further simplification is obtained assuming that 
the scale of variation of p{i',r) remains Oy = S~^^'^, as 
in the unperturbed case. We can then write p ~ po + pi, 
with 

51/2 

poiv, r) = fl-^poiv) = 27r^f^n ^^^^ 

the solution of Eq. (fT5|) with g{r/ry) set to zero, and 
estimate vdrPi ~ S^^/'^drPi, dy{vpi) — (l/2)i3o9^Pi ^ 
Vpi for some constant F. Substituting into Eq. 



[S~^'''dr + V]pi^ Bog{r/r,)pa. 



(19) 
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A new length S*"^/^, the distance travelled by a particle 
in a Stokes time, thus enters the problem. On the scale 
of S^'^/'^, the term g in Eq. p9|) behaves like a Dirac 
delta and we find: 



Pi 



{v,r) - ri-ie-i/2exp(-r5i/>|)poH. (20) 



Substituting into Eq. (|14p would lead to concentration 
fluctuations with 0(e^^/^) variance and 0{S^'^/'^) corre- 
lation length. 



V. THE VELOCITY STRUCTURE 

The analysis carried on in the previous section did not 
answer the question of which particle pairs arc responsi- 
ble for the production of the concentration fluctuations 
at large e. We are going to show here that concentration 
fluctuations are indeed produced by slowly approching 
particle pairs. 

We prove that the contrary is impossible and assume 
therefore that the dominant contribution is due to parti- 
cles approaching ballistically at scale r^, which spend a 
time T ~ e~^/^r5 at that scale. In order to show that the 
assumption of dominance of ballistic particles cannot be 
true, it is sufScient to calculate p{r) = / p{v,r)Av as the 
stationary solution of the evolution equation 

p{r- 0) = / dr'Av'pir- 0|i^', r'; -T)p{y' , r'; -T), (21) 



where p{v^ r; t) is a generic non-equilibrium PDF and 
p{r-Q\v',r'--T) = p{r{Q)=r\v{-T)^v' ,r{-T)=r') is a 
transition PDF for the dynamics described by Eqs. ^ 
|9]). If we take t <Si T <Si ts and \r\ < r^, for most 
particles [i.e. for typical values of vl—T)], we will have: 



|r(-r)| ^ |r(0)-K-r)T| »rv 



(22) 



which descends from the fact that, for T <^ ts and ~ 
ay, V can be approximated as constant. But Eq. ([22]) 
tells us that the dominant contribution to the integral in 
Eq. ([2T|) for p(r; 0) = p(r) (the equilibrium PDF), is from 
values of r' for which p(v' , r', — T) = p{v' ■, r') — po{v' , r'), 
where po is the spatially homogeneous PDF of Eq. ([18]). 
For the ballistic particle pairs (which we are assuming to 
dominate the statistics), we can calculate the correction 
to ballistic motion, writing r{t) — ro(t) + ri{t), where, 
using Eq. ([8]): 



ro{t) =r{-T) + u{-T){t + T), 
r^it) ^/!^rdr6(ro(r))e(r). 



(23) 



Using Eq. ([^5]) . we can express the transition PDF 
p{r{0)\v{~T),r{~T)) in Eq. ([21]) in terms of the 
equivalent PDF p{ri{0)\i^{~T),r{-T)) and, using the 
first of Eq. ((23]) . in terms of the conditional PDF 
/9(ri(0)|i^(-r),ro(0)). Substituting into Eq. ([HI, with 



p(r;0) = p{r) and p(i.',r';-T) ~ po{>^',r') = n-^po{j^'), 
we obtain: 



pir)^n-' / dri(0)di.(-T)po(^^(-T)) 



xpin{OM-T),rom 



roiO)=r-ri{0) 



(24) 



where use has been made of the relation p{r) = 
J dro(0)dri(0)p(ro(0),ri(0)),5(ro(0)-Hri(0)-r). We can 
Taylor expand Eq. (|24p around ro(0) = r and we obtain 



p{r) 



1 



where the average is on i'{^T) (in terms of po) ^md ^(r). 
Using Eqs. ([23]) and ([9]), we can write: 

p{r)^^dl{j\^ATB{r + uT)), 

where the remaining average is on v] setting r = 0, we 
see from Eq. ^ that the time integral is dominated by 
min(r, /\v\) and we obtain the final result 

p{r = 0) ^ n-\-^Bm{r./Hf\ W\ > r^/T), (25) 

where the conditional average to right hand side is peaked 
at ~ Tv/T and diverges for e^T/rg — > oo. 

Thus, ballistic pairs do not dominate the dynamics. A 
calculation taking into account small velocities in con- 
sistent way [S^, in turn, would lead to convergent ex- 
pressions for p(i/, r) and allows to show that the con- 
tribution to concentration fluctuation is concentrated at 
v ^ ^ (T^I^Uy, as could have been guessed from the 
analysis in Section III. 



VI. CONCLUSION 

The large Stokes number limit of inertial particles in 
a random velocity fleld is plagued by several subtleties. 
As pointed out in JJ, 20], when this limit is reached, the 
particle dynamics becomes dependent on a new parame- 
ter defined in terms of the Stokes and Kubo numbers S 
and K [see Eq. ([3])], namely: e — K'^S. This becomes ap- 
parent in the case of short-time correlated velocity fields, 
such as the Kraichnan model [s^, in which K = 0, and 
the role of the Stokes number is played indeed by e [from 
Eq. ([3]), is trivially equal to infinity]. Thus, as pointed 
out in [17.] , the large Stokes number asymptotics is really 
the large e one of Eq. ([TO]) . 

The point of this paper is that, for e 1, the concen- 
tration dynamics is not governed by particle pairs stay- 
ing close long enough for a local theory on the lines of 
[3,[2l],[3lj to work. While, for small e, the particle phase 
is locally mono-disperse in velocity [3g, for e > 1, the 
particle velocity distribution has finite width and only 
very few such almost co-moving particle pairs exist (see 
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discussion at the end of Section III). This property is 
expected to hold also in more than ID. 

Small relative velocities still appear to be important: 
as suggested in [28], the particles behave like a gas in 
thermodynamic equilibrium, but the correlations origi- 
nate, rather than from random collisions, fr'om particle 
encounters at small relative velocities, as made clear by 
the analysis leading to Eq. (I25|) . A qualitative argument 
in Section III, confirmed by the analysis in [s^], suggests 
that the only particle pairs that have a chance to be cor- 
related are those travelling at relative velocity ^ e~^/®tT„, 
that is the minimum for ballistic relative motion at scale 
Tv. However, these velocities are still 0(exp(ei/3)) lar ger 
than those producing the long time dynamics needed for 
clustering in the e ^ 1 limit. 

The concentration fluctuation variance appears to de- 
cay rather slowly, like e~^/^ at large e [see Eq. pO|]. This 
decay rate could be interpreted as the product of the con- 
centration fraction e~i/6 of the particles staying close 
long enough for their relative velocity to be modified in 
some way, and their permanence time ^ g^i/^ (in units 
Ts) at separation |r| < [see Eq. p^ ]. 

This mechanism of concentration fluctuation produc- 
tion leads to an important modification of the relative 
velocity distribution. As particle clusters are the result 
of particle pairs with small relative velocities, the total 
velocity PDF will be weighed more to small relative ve- 
locities. This results in a depletion of the relative veloc- 
ity variance, quantified in the exact relation (|17p . that 
becomes substantial at e 1. It is to be ascertained 
whether this property is preserved in more than ID. 

Summarizing, we have the following situation: extrap- 
olating the small e theory of [2l|,[2l,[23| to e > 1, predicts 
that clusters become less and less singular, making even- 
tually a transition from fractal to space-filling objects. 
Extrapolating the opposite way around, the present large 
e theory predicts "residual" (not singular) concentration 
fluctuations, becoming of the order of the mean concen- 
tration as e — > 1, due to more and more particle pairs 
contributing to the fluctuations. It is reasonable to ex- 
pect that the situation at e 1 will be a combination of 
the two pictures. 

In dill, it was suggested that clustering is not cen- 
tral in determining collision rates of relevance e.g. in rain 
formation, as instead are caustics formation and a finite 
width of the velocity distribution. If the present picture 
continues to hold in more than ID, it would result in 
an even stronger statement, that is, at large e, cluster- 
ing and modifications to the velocity distribution work 
against one another. This is in some way analogous to 
the situation at e — > 0, of clustering dominance accompa- 
nied by caustics disappearance and vanishing of collision 
rates [2l| (see also [37,]). 

Another feature of concentration fluctuations at large 
e, is the strong scale separation between their decay 
length and the correlation length of the fluid: the first is 
by a factor ~ e^^^ larger than the second [see Eq. (|^D|) ]. 
The concentration fluctuation decay length is basically 



the distance travelled by a particle in a Stokes time and 
is independent of the fluid correlation structure. Taking 
into account also the e^^/^ variance decay, smaller vor- 
tices may thus contribute in a non negligible way to the 
concentration fluctuations of particles with Stokes time 
in the turbulent inertial range. 



APPENDIX A: FINITE STOKES NUMBER 
CORRECTIONS 

The Langevin approximation to Eq. ^ provided by 
([5]) is the lowest order in an expansion in inverse pow- 
ers of the Stokes number S. Let us calculate the next 
order. Given initial conditions at i = 0, Eq. ^ can be 
integrated to give: 



v{t) = v{0)e- 



dt'e''-'u{x{t'),t') 



and we can write x{t) = x(t) + x(t) with 



x{t) = x{0) + v{0){l 



and 



xit) 



dt'{l 



u{x{t'),t') 



Following an approach similar to 6], the perturbation 
expansion is carried on with respect to x and the first 
correction to v(t) is given by, for t <ti ts = 



v{t) = t;(0)(l -t)+ /g dt'u{x{t'), t') 

+ /o dt' /„*' dt"{t' - t")u'{x{t'), t')u{x{t"),t") 



{Al) 



where u'[x, t) — dxu[x, t) and x{t) ~ x[Q)+v{Q)t. Taking 
t :S> te, we find that the second piece to right hand side 
of Eq. (Al) behaves like a Wiener increment; using Eq. 



dt'uix{t'),t') = 2Tp{v{0))t 



and we recover Eq. ([5]). Notice that the substitution 
x{t') — > x{t') avoids a definition of the noise in implicit 
form and ^(t) in Eq. ([5]) is automatically defined in the 
Ito prescription [sH . 

The last term in Eq. (Al) contains a drift correction 
in the form, from Eq. ([1]): 

{u'{x{t'),t')u{x{t"),t")) = -{KS)MO)it' -t") 
x{u{v{0){t' ^t"),t' ~t")u{0,0)). 

Equation (Al) takes then the final form, in the limit t 
dt: 

dv = -[1 + -f{v)]vdt + {2Tp{v)y^^dw {A2) 



where 



7(w) = {KSf / dt t^{u{vt,t)u{Q,0)). 
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and w{t) = Jp ^(T)dr is the Wiener increment. The noise 
amphtude is defined up to 0(dt^/^) terms associated with 
the fluctuating part of the second line of Eq. (Al). 

Thus, the first order correction to the Langevin equa- 
tion for v{t) is an 0{K^/S) renormahzation of the Stokes 
time, whose sign depends on the profile of the correla- 
tion function for u{x{t),t): a strictly positive correlation 
corresponds to a decrease of the Stokes time, while oscil- 
lations in the correlation may lead to an increase. This 
would lead to a decrease (an increase) of the drift induced 
on the particle by an external force. 

Notice that the drift correction 7 provided in Eq. (A3) 
does not necessarily coincide with what would be ex- 



pected from a change from Stratonovich to Ito prescrip- 
tion [3I]. For instance, adopting for the random field 
correlation the explicit expression g{x) = exp(— a;'^/2), 
would lead to twice what would be obtained interpret- 
ing the substitution x{t') — > x{t') in dt'u{x(t'),t') as 
a Stratonovich to Ito prescription change. Indicating 
by ds the Stratonovich differential and writing b{v) = 
{2T,{v)f/^: 

h{v)(isw — —b{v)b' {v)dt + b{v)dw, 

we find in fact, from Eq. ([6]): ^b{v)b'{v) = —\"iv. 
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